Finite-size corrections of the Entanglement Entropy of critical quantum chains 
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Using the density matrix renormalization group, we calculated the finite-size corrections of the 
entanglement Q-Renyi entropy of a single interval for several critical quantum chains. We considered 
models with (7(1) symmetry like the spin-1/2 XXZ and spin-1 Fateev-Zamolodchikov models, as well 
as models with discrete symmetries such as the Ising, the Blume-Capel, and the three-state Potts 
models. These corrections contain physically relevant information. Their amplitudes, that depend 
on the value of a, are related to the dimensions of operators in the conformal field theory governing 
the long-distance correlations of the critical quantum chains. The obtained results together with 
earlier exact and numerical ones allow us to formulate some general conjectures about the operator 
responsible for the leading finite-size correction of the a-Renyi entropies. We conjecture that the 
exponent of the leading finite-size correction of the a-Renyi entropies is ~ 2X^1 a for a > 1 and 
pi = V, where denotes the dimensions of the energy operator of the model and v = 2 for all the 
models. 

PACS numbers: 03.67.Mn, 05.50.-(-q, 05.65. -|-b 



I. INTRODUCTION 

Usually, we are interested in ground state properties 
of quantum critical chains in the thermodynamic limit. 
However, obtaining the physical properties in this limit, 
in general, is not possible. On the other hand, many nu- 
merical techniques that give accurate information about 
the system, work for finite systems. It is therefore of fun- 
damental importance to obtain the physical properties of 
the infinite system of interest from its formulation in a 
finite-size geometry. Indeed, various physical quantities 
of interest are directly related to the finite-size scaling 
corrections of some physical measure. In this work, we 
are interested in the scaling corrections of the entangle- 
ment entropy of one-dimensional critical systems, which 
has been intensely debated in recent years. 

Consider a quantum chain with L sites, described by 
a pure state whose density operator is p. Let us consider 
that the system is composed by the subsystems A with 
I sites {I — 1, • ■ • , and B with L — I sites. The Renyi 
entropy is defined as 



The first term in this equation, is the conformal field 
theory (CFT) prediction and is given \>y^— 
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where c is the central charge and da is a non-universal 
constant. 

Laflorencie et al in Ref. @were the first to notice, in an 
investigation of the spin-1/2 Heisenberg chain under open 
boundary condition, that unusual corrections to scaling 
appear in the Renyi entropy Sa{L,£) with a = 1. They 
noted that the standard CFT term [Eq. |3] could not ex- 
plain a strong oscillation observed in the von Neumann 
entropy, which seems to have origin its in the oscillating 
behavior of the spin-spin correlations of the Hamiltonian. 
These strong oscillations, for the quantum open chains, 
were observed subsequently by several authors.—"— More 
recently, Calabrese et al. in Ref. |^ (see also Ref. [13) in- 
vestigated the anisotropic spin-1/2 Heisenberg chain with 
PBC and they noted that although those oscillations are 
not present in the von Neumann entropy Si, they are 
still present for the Renyi index a >1. Based on exact 
results and also in numerical calculations of the spin- 
1/2 XXZ chain with PBC and zero magnetic filed, Cal- 



where pj^ = Tr^p is the reduced density matrix of the 
subsystem A. The von Neumann entropy, also known as 
entanglement entropy, is the particular limit a ^ 1. 

In the scaling regime 1<<£<<L, itis expected that 
for the critical one-dimensional systems, under periodic 
boundary conditions (PBC), the Renyi entropy of the 
ground state behaves as 



abrese and collaborators conjectured that S, 
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from a non- universal constant get has the following uni- 
versal behavior§rii(see also Ref. Q) 
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where for convenience we introduce the function 
Da{L,£). In (4) Pa is a critical exponent and k, which 
has distint values for different models, gives the spatial 
period A = — of the oscillations Jj'^ For instance, for the 
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spin-s XX Z chains at zero magnetic field kxxz = t^- 
For these chains, the universal exponent is Pa = 
where K is the Luttinger liquid parameter of the under- 
lying CFT.~iii- For the Ising model nising — 2kxxz, and 
we see no oscillations. The origin of the oscillating factor 
[cos(k£)] was not yet completely understood, however it 
has been observed in the systems whose spin correlations 
show an oscillating behavior. 

By now, the observation of the unusual corrections to 
scaling in entanglement entropy with the predicted expo- 
nents [Eq. were confirmed only for: the spin-s XXZ 
chainS)2iii the ID attractive Hubbard model, and in a 
dipolar boson quantum chainJ^ As observed by Cardy 
and Calabrese in Ref. [l^ the origin of the exponent Pa 
in Eq. ^ is the conical singularities produced by the 
conformal mapping used to describe the reduced density 
matrix — Ti^p in CFT. This exponent is related to 
the scaling dimension X'^"" of an operator of the under- 
lying CFT by Pa = 2X''°"/aM. This operator can be 
relevant {X"""^ < 2) or not {X"'"' > 2). Beyond these 
corrections, we may also have the usual ones coming from 
the irrelevant operators (dimension X^ > 2), whose lead- 
ing contributions to ^ arc of order for any 
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Recently, Calabrese and Essler in Ref. [lO], based on 
an exact calculation for the XX quantum chain, found 
the appearance of a correction term of order £~'^ for any 
value of a. As discussed in Ref. [H, we should expect 
a-independent correction terms from the irrelevant op- 
erators. This would imply that the contributing irrel- 
evant operator has dimension X^ ~ 3. Although such 
an irrelevant operator exists in the conformal tower of 
the XX model, it is quite unlikely that this contribution 
comes from such an operator. In this model, the oper- 
ator responsible for the leading finite-size corrections of 
the eigenenergies has dimension X^ = which would 
produce a correction term of order £~'^. A possible expla- 
nation is that the conical singularities may also produce 
an a-independent contribution to the corrections to scal- 
ing, besides the term of order ^-2X'=°"/q_ Although not 
completely clear, such contribution might come from the 
combination of contributions of the relevant operators^^ 
and in the case of the XX chain it ends up with the main 
contribution of order £~'^ , with ly = 2. 

Although the CFT analysis predicts the possible finite- 
size corrections of the Renyi entropies, several general 
questions remain to be answered for a better understand- 
ing: a) Among the possible relevant operators in the 
CFT what should be the leading one, with dimension 
responsible for the corrections coming from the 
conical singularities giving a-dependent contributions of 
order ^-2-^"°"/" ? What should be the leading irrel- 
evant operator with dimension X^ , giving the contribu- 
tion 0{i~^'-^ ^^^) ? c) Are the a-independent correc- 
tions £~'^ , coming from the conical singularities, as hap- 
pened in the XX case, general ? In this case what is the 
value of ? d) The oscillating behavior observed in the 
a > 1 Renyi entropy is absent in the quantum Ising chain. 



Should we obtain the oscillatory behavior of the entropies 
for other quantum chains with C/(l) symmetries, like the 
XXZ chain? 

In order to test the above predictions and answer the 
above questions, we present in this paper the calculation 
of the corrections to scaling in the entropies of several 
critical quantum chains belonging to distinct universal- 
ity classes of critical behavior. We consider models with 
Z{2) and Z(3) symmetries. The models with U(l) 
symmetry are the spin- 1/2 XXZ quantum chain and the 
spin-1 Fateev-Zamolodchikov model. — The models with 
Z{2) and Z{3) symmetries we consider are the quantum 
Blume-Capel and the quantum 3-state Potts model, re- 
spectively. 



II. RESULTS 

From the discussions presented in the introduction we 
expect, in the region where the system size L and subsys- 
tem size i are large (L >> £ >> 1), the general behavior: 



Da{L,£) 









sm — 


+ ba 


Ht)] 







-2(X'-2) 



■(5) 



The second and third terms are the leading contributions 
due to the conical singularities and the last term is the 
leading contribution due to the standard correction to 
scaling operator, with dimension X^ . The operator that 
produces the second term in (6) has dimension AT"^"", 
while the relation of the a-independent exponent v with 
the dimension of operators is unknown. Models exhibit- 
ing oscillations in the a-entropies have 7^ and k ^ 
(mod 27r). Note that only the exponent of the second 
leading corrections depend on the value of a. We intend 
to evaluate the power of the dominant term in ([5]) , which 
we denote by Pa ■ In order to estimate the exponent Pa , 
we fit our data with the following equation 



Da{L,£) = da+fa[cos{K£)f *° 
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We report first the results already known from earlier 
studies for some models. The spin-1/2 XXZ chain shows 
oscillations for a > 1. These corrections come from the 
operator with dimension X'^"" = K , where K is the Lut- 
tinger liquid parameter.^ This is the scaling dimension of 
the energy operator of the modelj^ For a = 1 the oscilla- 
tions do not appear, i. e., 171 = 0, and the dominant term 
in ([5]) is not known for general values of the anisotropy of 
the model. In the particular case where we have the XX 
model (free fermion case), the dominant correction for 
Di{L,£) is given by the third term in (6) with ly = 2.^" 

The non integrable spin-s XXZ chain (s — 1, 3/2, . . .) 
on its critical regions also shows entropy oscillations for 
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a > 1, but the amplitudes ga decreases strongly as 
we increase the value of the spin sJ^ On these mod- 
els, the amplitudes of the entropy oscillations are also 
ruled by the energy operator of the model, with dimen- 
sion = = K, given by the Luttinger parameter 
of the underlying CFT. 

The results for the Ising model can be obtained from 
those of the spin- 1/2 XX model due to the exact corre- 
spondence of their entropies, as shown in Ref. [l!]. The 
Ising model does not show oscillations, for any a, due 
to the fact that nising = '^i^xxz = 27r. We also have 
gi — Q and (^q, ^ (a > 1) with the dominant correc- 
tion term given by = 1, which is also the dimension 
of the energy operator of the model. In the case a = 1 
the dominant correction comes from the third term of (5) 
with u = 2. 

Below, we present new numerical results for three mod- 
els: (i) the Blume-Capel Model (BCM), (ii) three-state 
Potts Model (3SPM), and (iii) the Fateev-Zamolodchikov 
quantum chain (FZQC). We investigated these models 
with the density-matrix renormalization group (DMRG) 
technique^ under PBC. We have done 4-8 sweeps. For 
the BCM and the 3SPM we used typically m = 400 states 
per block. This number of states kept in the truncation 
process is enough to give very precise results, the dis- 
carded weight being typically about 10""'^^. However, for 
the FZQC we used a much larger number of states in 
order to obtain precise results (up to m = 3000). In this 
case, the discarded weight was typically 10 in the final 
sweep. 



00 



(a) 



0.364 



0.357 
10 



'IM a=l A.=1.0 
•IMa=3 A.=1.0 
.BCMa=l 1^1.1 
•BCM a=3 1^1.1 




L=96 



20 



40 



(b) L=96 








' /• ISING 


?i=1.0 a=3 


/ " BCM 


7=1.1 a=3 







20 



30 



40 



Figure 1: (Color online). Results of the Renyi entropy 
Sa (L, £) and of the function Da {L, I) for the Ising and Blume- 
Capel models with PBC and size L = 96. (a) Sa{L,£) vs £ 
for two values of q, and some values of coupling (see legend). 
The symbols are the numerical data and the solid lines con- 
nect the fitted points using Eq. (b) Da{L,£) vs £ for 
Of = 3. We added 0.03 in the values of the Ising entropy in 
order to see better both data in the same figure. 



A. The spin-1 Blume-Capel model 

We begin by introducing the spin-1 BCM quantum 
chain. This model is obtained by the time-continuum 
limit of the well known BCM in two dimensions.— It de- 
scribes the dynamics of spin-1 localized particles, with 
Hamiltonian given by 

i 

where and are the spin-1 SU{2) operators. The 
phase diagram, in the 5 — 7 plane, is known from earlier 
numerical finite-size scaling studies based in the crossing 
of the mass gap energies (see Fig. 1 of Ref. [2^. For 
values 7 > 7tr the Hamiltonian has a quantum critical 
line Sc{'y) governed by a CFT in the same universality 
class of the quantum Ising chain, i. e., central charge 
c = 1/2. At 7tr the model has a quantum tricritical 
point at StT in the universality class of the tricritical Ising 
model, having central charge c = 7/10. For 7 < 7tr 
there is a line S = Sgap{j) of first-order phase transitions. 
Recently, highly accurate estimates of some points of the 
critical line were obtained by a new method based on 
the entanglement entropyi^ In particular, the tricritical 
point was located at j^j- = 0.41563 and S^j. = 0.91024. 



Since the model has a critical line in the universality class 
of the Ising model (IM) it is interesting, for the sake of 
comparison and as benchmark tests, to consider also the 
integrable IM. The Hamiltonian of the IM is given by 

j 

where tr^, cr^ are spin-1/2 Pauli matrices. This model 
has a critical point Xc — I that can be obtained from 
its exact solution, or even more simply from its self dual 
propertyjSSi 

In Fig. [IKa), we show the Renyi entropy Sa{L, £), a,s a 
function of i, for the IM and the BCM with PBC and lat- 
tice size L = 96. In the IM case, we select the exact criti- 
cal point A = 1. Since in the BCM we do not know the ex- 
act values of the critical points, we used the best estimate 
of the critical points available from previous work.-^ In 
particular, for 7 = 1.1 we used Sc = 0.31357, where we ex- 
pect an error in the last digit (see Ref. 24). In both cases, 
the fits to Eq. ^ give a central charge value very close to 
the exact one. We found that jc/** -c'=^°'=*| < IQ-^ for aU 
critical couplings considered (some of these estimates are 
presented in Table |TT| . We emphasize that the discrep- 
ancies between the exact values and the numerical data 
are due to the finite-size of the systems considered. The 
numerical errors (inferred from the dependence of the en- 
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tropy with m) are smaller than the size of the symbols 
in the figures. For the models with discrete (continuous) 
symmetry the errors are smaller than 10~^ (10^^). 

The results presented in Fig. [TJa) are nice illustra- 
tions of the universal behavior of the Renyi entropies at 
the critical points. Although the entropies are calculated 
from ground states defined in Hilbert spaces of quite dis- 
tinct dimensions (2^ and 3^) the entropies show the same 
behavior. The small deviation shown in this figure is due 
to the non- universal constant da in (|3]). Notice also that 
both models, in the Ising Universality class, show no os- 
cillations in the entropies (even in the OBC case, not 
shown), i. e., K = in ([S]). 

The corrections to scaling contributions are calculated 
from the behavior of Da{L,i) given in ([S]). As illustra- 
tive examples, let us consider, as in Fig. 1(a), the IM 
at the critical coupling = 1 and the BCM at the cou- 
plings 7 = 1.1 and Sc — 0.31357. In Fig. [Hb), we present 
Da{L, i) for the critical couplings of the quantum chains 
with PBC, a = 3, and lattice size L = 96 (similar results 
are found for L = 64). Since at these couplings, both 
models have central charge c = 1/2 we set this value in 
Eq. ([3]). Note in these two figures that the changes in 
Da{L,£), as function of £, are quite small, which indi- 
cate that the non- universal constants ga , Oq, and ba that 
appear in ([SJ are very small. This fact makes the task 
of determining the exponent in ([6]) a challenge. Since 
we are interested in the asymptotic behavior of Da{L, £), 
we have discarded the first 10 values of £. Comparison 
of the entropies with different values of m indicate that 
the errors in the DMRG evaluations of the entropies are 
around ~ 10^^ — 10^^. Due to this fact, we also discarded 
the sites £ at which \Da{L,£) - Da{L,£+l)\ < 10"^. 

We get Pa by fitting our data to Eq. ([6]). The solid 
lines, in Fig. [T](b), connect the fitted points. In order 
to check if we are able to extract good estimates of pa 
with the above mentioned procedure, we compare the 
values of Pa-, acquired with the fit procedure with the 
exact ones {pa = ^) for the Ising model. The finite-size 
estimates for the exponents P2, Ps and p^ are depicted 
in Table ID As we can observe in this table, the finite- 
size estimates differ very little (less than 5 x 10^'^) from 
the exacts values. We then believe that this approach 
will also provide good estimates of the exponents pa for 
other models like the BCM. 

Before presenting our results for the BCM, it is inter- 
esting to mention that, for the IM, we got only 8 points of 
DalT{L,£) satisfying \Di{L,£) ~Di{L,i+l)\< 10"^ 
for £ > 2. This in fact is expected, since for a = 1 in 
the IM, the dominant correction to scaling is conjectured 

as Di{L,£) =di + io [LsinC/L}V + "^^^^^ ^ 

larger exponent than the ones of a >1. It is amazing that 
a fit of those 8 points to Eq. ^ give the following values: 
di = 0.478, /i X (L/nf = 0.018 and pi = 2.05. These 
values are very close to the exact ones, i. e.: di = 0.4785, 
/i X (L/7r)2 1/60 0.0166 and pi = 

It is interesting to stress that the correction term of 
order £~'^ in the case of the IM does not come from the 



P2 P3 P4 

0.978 0.646 0.473 
(1.000) (0.6666) (0.5000) 



Table I: The finite-size estimates of the exponents obtained 
from the fit of Eq. ^ for the Ising model with PBC and 
L = 96. The results in parentheses are the exact values. 



second term in ([5]) with Kjsing = 27r, but from the third 
term with v = 2. This is due to the exact relations among 
the entropies of the IM and the XX chain, and the fact 
that there is no oscillations in the last model for a = 1, 
i. e. 51 = 0. 

Let us now focus on the BCM. We show in the last 
column of Table HIl the finite-size estimates of the central 
charge c for three values of 7 along the critical line. We 
got these estimates by a simple fit of our data to Eq. 
^ for a = 1 . As we note, theses values of 7 satisfy 
7 > 7tr = 0.41563, and the system belongs to the same 
universality class of critical behavior as the Ising model 
with c = 1/2. Following the same procedure done in the 
IM, we fit our data to Eq. ^ with the appropriated val- 
ues of c [see Fig. [^b)], in order to extract the exponents 
Pa- Some finite-size estimates of are depicted in Ta- 
ble [TTl As observed in this table, for 7 > 7tr, our results 
are consistent with the exponent Pa = 2 /a, which is the 
same of the IM. In the case a = 1, where the evaluations 
are more difficult we got, in the BCM, around 10 points 
of Da(L,e) satisfying \Di(L,£) - Di(L,£+l)\ < 10"^ 
for ^ > 5, and the fit to Eq. ([6]) indicate that pi — 2. 
Since the energy operator — I, this term could be 
produced either by the first term in ([5]) with gi ^ but 
K = 0, or by the second term in ([5]) with ly = 2. However, 
the exact result derived for the IM, that we expect to be 
the same for any model on its universality class, indicate 
that gi — and the dominant term in ^ is the one with 
exponent v = 2. Our estimates for this exponent pi = 1/ 
are shown in Table HTl 

At the tricritical point ('jtric, Stric) ^ we were not able 
to extract the exponents Pa due to the fact that the a- 
entropies are very sensitive to the coupling constants of 
the model. We observed that a very small error (e ~ 
10~^) in the critical coupling Sc affects very little the 
value of the entropy along the critical line. Typically, 
the entropy changes are of the same magnitude of e. On 
the other hand, at the tricritical point, a small change in 
the critical couplings produces a much larger effect. Due 
to this fact, since the precision of the tricritical couplings 
are ^ (10~^), which is of the same order of magnitude as 
Da, we are not able to obtain a reasonable estimate of 
Pa at the tricritical point. 
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_7 P2 P3 Pi V e 

1.2 1.066 0.697 0.517 2.071 0.5000 

1.1 1.041 0.682 0.504 2.062 0.5001 

0.7 1.154 0.776 0.589 1.955 0.5015 

(1) (2/3) (1/2) (2) (1/2) 



Table II: The finite-size estimates of the exponents pa ac- 
quired by the fit of Eq. (O for the BCM with PBC, L = 96, 
and some values of 7. In the last two columns we also present 
the finite-size estimates of the conformal anomaly c, and the 
exponent p\ = v (see text). The values in parentheses are the 
expected ones. 



B. The three state Potts Model 

The quantum 3SPM is a quantum chain obtained by 
the T-continuum limit of the two dimensional classical 
3SPM.^* The quantum Hamiltonian of the model can be 
written as 

Hpotts = - ^(i?^i?l+i + h.c.) + AO,, (9) 

i 

where R and O are 3x3 matrices given by 



"3Sp1vi a=l 
"3SPM a=3 
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Figure 2: (Color online). Renyi entropy Sa[L,t) and the 
difference (L, i) for the 3SPM with lattice size L = 96. The 
symbols are the numerical data and the solid lines connect the 
fitted points (see text), (a) Sa{L,i) vs £ for A = 1, and some 
values of a (see legend), (b) Da{L,(.) vs £ for A = 1, a = 2 
and Q = 3. 
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As the IM, due to its self duality the 3SPM also has 
the critical point at Ac = 1. However, the critical fluctua- 
tions of the 3SPM are described by a CFT with c = 4/5. 
For PBC, as a consequence of the modular invariance of 
the related two dimensional model defined on a cylin- 
der, the associated CFT is described in terms of 10 pri- 
mary operatorsi^ Among these operators, we focus on 
the most relevant operators, which may lead to the un- 
usual corrections to scaling in the entropy. 

The first nonzero dimension (X^ = 2/15) is associated 
to the order parameter and the second one {X^ = 4/5) is 
associated to the energy operator. The lowest irrelevant 
operator responsible for the finite-size corrections of the 
eigenenergies of the Hamiltonian has dimension — 
14/5^ 

Note that in the IM and the BCM, the leading correc- 
tions to scaling in the a-entropies with a > 1 are related 
to the dimension of the energy operator . For this rea- 
son, we expected in the 3SPM that pa = = for 
a > 1. Below we use the same procedure used in the IM 
and the BCM to get the exponent Pa- 

In Fig. [2ja) , we present the Renyi entropy as a function 
of £ for the 3SPM at the critical coupling Ac = 1 for a 
system with lattice size L = 96. As observed in this 
figure, we also do not see the parity oscillations, like the 



models considered in the last section (which also have a 
discrete symmetry). We are able, in this case, to obtain 
a nice fit of our data by considering only the standard 
term of the CFT [Eq. ([3])]. The solid lines in these figures 
connect the fitted points. In this case, we get c — 0.799 
(see Table Hill) very close to the expected value c = 4/5. 

Finally, we depicted in Table IIIII the exponents pa ob- 
tained through a fit of our data to Eq. As we can 
see in this table our numerical results indicate that, for 
a > 1, Pa = confirming our prediction that like in the 
IM and BCM the energy operator gives the most impor- 
tant contribution in the conical singularities. In the case 
of the von Neumann entropy, i. e. a = 1, our numeri- 
cal results indicate that the leading corrections are given 
by the third term in (6) with the parameter v — 2. As 
happened in the Ising universality class the operator that 
gives the leading contribution for a > 1, that we believe 
to be the energy operator, does not contribute for a = 1, 
i. e., gi = 0. Moreover, since the leading irrelevant op- 
erator is X^ = 14/15, we would expect a contribution of 
order [sin(7r^/L)]~^^^^ in which is also not present, 
i. e., 61 = 0. This implies that both the energy and the 
leading irrelevant operator ruling finite-size corrections 
do not contribute, at least at leading order, to the usual 
von Neumann entropy. The result gi — found for all 
models discussed above indicate that, probably, this is a 
general behavior for quantum critical chains. 



6 



P2 



P3 



Pi 



0.829 0.523 0.371 0.799 
(0.800) (0.5333) (0.400) (0.800) 



Table 111: The finite-size estimates of the exponents Pa ac- 
quired by fit of Eq. ((6]) for the three state Potts model for 
L — 96. In the last column we also present the finite-size 
estimates of the conformal anomaly c (see text). The results 
in parentheses are the expected values. 
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C. The spin-1 Fateev-Zamolodchikov quantum 
chain 
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The spin-1 Fateev-Zamolodchikov quantum chain 
(FZQC) is an exact integrable quantum chain whose 
Hamiltonian is given hy^ 


CO 
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7 ■ 
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a=2 PZQC 
(x=3 



Hpz — ( 



j 



-4sin2(V2) {T^T^ + h.. 



rj-iZ 



(10) 



where e 



±1, Tf 



and 



s|, s^, s| are the spin-1 SU{2) operators. 

This Hamiltonian has a J7(l) symmetry, having the z- 
magnetization as a good quantum number. The antifer- 
romagnetic (e = 1) and ferromagnetic (e = —1) models 
show a critical phase with continuously varying expo- 
nents for < 6 < 7r/2. The antiferromagnetic and fer- 
romagnetic models are governed by a CFT with central 
charge c = 3/2 and c = 1, respectively.^"^ The CFT 
in the c = 3/2 case is described in terms of composite 
fields formed by Gaussian and Ising operatorsi^— The 
energy operator, which we believe is responsible for the 
leading contributions due to the conical singularities, has 
dimension = ^^^2S) ~^ I leading irrelevant 

operator responsible for the finite-size corrections of the 
eigenenergies has dimension = ^^^^ + li^i^ 

In the case c = 1, we have a Gaussian CFT or a 
standard Luttinger liquid, whose Luttinger parameter 
changes continuously along the critical line. In this case, 
the dimension of the energy operator is = ^ and the 
leading irrelevant operator responsible for the corrections 
to scaling has dimension X^ = "^-^^ 

In Figs. [3l[a)-(b), we present the Renyi entropy 
Sa{L,t} for spin-1 FZQC with, PBC, 5 = 0.5, and size 
L = 72. Notice that for e = 1, differently from the Ising 
and the Blume-Capel models, the a— Renyi entropies os- 
cillate for a > 1, which is a signature of the oscillating 
spin-spin correlations of the quantum chain. The absence 
of oscillations for a = 1 also suggests that like the other 
models we have studied in previous subsections 51 = in 
Let us focus, first, on the case a = 1. In this case, we 
get very nice fit of our numerical data (black triangles) to 




Figure 3: (Color online). Results of the Renyi entropy 
Sa{L,£) and the difference Da{L,£) for the spin-1 Fateev- 
Zamolodchikov quantum chain with PBC and size L = 72. 
The symbols are the numerical data and the solid lines con- 
nect the fitted points (see text), (a) Sa{L,l) vs I for e = 1, 
5 = 0.5, and some values of a (see legend), (b) Same as (a) 
but for e = —1. (c) Da{L,£) vs £ for e = 1, a = 2 and q = 3 
at (5 = 0.5. 



the standard conformal field term S"^^^ [Eq. The 
fit gives us c = 1.491 (c = 1.02) for e = 1 (e = -1), which 
is very close to the exact value Cexact = 3/2 {cexact = !)■ 
Similar estimates of c were also obtained by considering 
a = 2 or a = 3. Note that the central charge c of the 
FZQC is larger than the ones of Ising and Blume-Capel 
models. This means that the subsystem A is more en- 
tangled with the subsystem B in the FZQC, as compared 
with the critical Ising or Blume-Capel models. Due to 
this fact, we have to keep a much larger number of states 
TO (typically around m ~ 2000 - 3000) for the FZQC, in 
order to get results with some similar accuracy. 

Now let us consider the antiferromagnetic model e = 1 
and a > 1. In this case, it is not possible to fit the 
numerical data if we consider only the standard CFT 
term S^^'^ [Eq. ([3])]. Relative nice fits are obtained 
only if we consider besides the CFT term, the oscillatory 
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5 I P2 P3 Pi 

0.25 0.413 0.286 0.215 

(0.4223) (0.2815) (0.2112) 

0.5 0.432 0.304 0.231 

(0.4917) (0.3278) (0.2459) 

1.0 0.778 0.541 0.402 

(0.8129) (0.5419) (0.4064) 



Table IV: The finite-size estimates of the exponents pa ac- 
quired by fit of Eq. ((5]) for the spin-1 Fateev-Zamolodchikov 
quantum chain with PBC, L = 72 and some values of 5. The 
values in parentheses are the expected ones (see text). 



correction term [Eq. (U)] with n = 7r/2. The solid lines 
in Fig. are the fits to Eq. For instance, for 

a = 2 (a = 3) we get c = 1.48 (c = 1.45) and p2 = 0.54 
b3 = 0.36). 

Since we are confident that c = 3/2 for the spin-1 
FZQC, we set c = 3/2 in Eq. ©, as weh n = it/2 
in the function Da(L,£) [Eq. ([S])], in order to evaluate 
the exponent Pa in the same way we did in the Ising, 
Blume-Capel and three state Potts models. In Fig. [2Ib), 
we show Da{L,i), as function of i, for the spin-1 FZQC 
with PBC, e — 1, 6 — 0.5, L — 72 and three values 
of a. Even keeping m = 3000 states in the DMRG, it 
is very hard to get high precision in the values of the 
entropy for the spin-1 FZQC, as we already mention be- 
fore. Although the truncation errors are around ^ 10~^, 
comparison of the entropy results with different values 
of m, indicate that the errors in the entropy are around 
~ 10~^. For this reason, we discarded the points that 
\Da,{L,i) - Da(L,e + 2)1 < 10-^, as weh as the first 
10 sites. The solid lines in this figure connect the fitted 
points. The finite-size estimates of Pa in ([B]), obtained 
by this fitting procedure, are depicted in Table IIVI The 
results presented in this table indicate that Pa = 2X^/a, 
where = ^^Z.2S) ~^ ^ dimension of the energy 

operator of the model. This means that, like the other 
quantum chains, the energy operator gives the most im- 
portant contribution in the conical singularities for a > 1 . 
In the case where a = 1, like the other models, this con- 
tribution seems to be absent, i. e., gi — [see Fig. [2Ua)]. 
In this case, we were not able to estimate numerically the 
power of the leading correction, since we only have three 
points satisfying \Di(L,£) - Di(L,£ + 2)\ < lO'^. 

In the ferromagnetic case (e = —1) although we ex- 
pect oscillations to be present in the a-entropies with 
a > 1, we could not see them numerically [see Fig.[3l^b)]. 
This is consistent with the conjecture that the oscilla- 
tions are ruled by the energy operator. In this case 
the amplitude of the oscillations will decay like i'P", 
Pa = = The region where we would better 

see the oscillations would be for i5 « ^, where Pa has 



the smaller values. However the sound velocity, which 
IS given hy^ Vs = 7rsin(2(5)/(27r - 26) is close to zero 
in this region. This makes the convergence quite slow 
in the DMRG, due to the small size of the energy mass 
gaps. To avoid this problem and get enough precision, 
we considered values of d in other regions, like 6 = 1/2 
for example [see Fig. [3l^b)]. However, in this case, the 
expected value of Pa = 27r/a is for a = 3, Pq, w 2.1. 
This gives a strong decay, large enough to forbid, within 
the numerical precision we have, the observation of the 
oscillations for ^ > 10. 



III. CONCLUSIONS 

Most of the critical quantum chains are conformal in- 
variant. This symmetry implies that the mass gap ampli- 
tudes of these critical chains, in a finite lattice, are related 
to the conformal anomaly and critical exponents that la- 
bel the particular universality class of critical behavior 
Due to these relations, the most frequently used method 
to extract the conformal anomaly and critical exponents 
comes from the finite-size scaling of the eigenenergies of 
the critical quantum chains. 

On the other hand, the conformal anomaly can also be 
calculated from the Renyi entanglement entropies of the 
ground state of finite quantum chains ^i^i^i^^ a property 
related to the eigenf unctions instead of the eigenenergies. 
These results raised the question about the possibility to 
extract all the critical exponents from a complete finite- 
size study of the Renyi entropies of the critical quan- 
tum chains. This question has also practical implications, 
since in DMRG calculations the entanglement entropies 
are much simpler to calculate, as compared with the mass 
gap energies. Unfortunately the leading terms of the en- 
tanglement Renyi entropies of the low-lying excited states 
are the same as that of the ground statei^i^ The crit- 
ical exponents only appear in the finite-size corrections 
of these entropies.^' ^® In the finite-size study of the en- 
ergy gaps the exponents are given by the leading terms 
and we know precisely the correspondence between the 
mass gap amplitudes and the dimensions of the under- 
lying CFT. In the finite-size study of the entanglement 
entropy, since the exponents are given not by the lead- 
ing term, but by the finite-size corrections, we do not 
have such exact correspondence. It is known from con- 
formal invariancei^ that these corrections are ruled either 
by relevant operators, due to the conical singularities in 
the conformal mapping, or by irrelevant operators. How- 
ever, it is not known what are those operators, or equiva- 
lently what should be the critical exponent that appears 
in the leading finite-size corrections of the entanglement 
entropies in general models. The results presented in 
this paper, together with known results induce us to an- 
nounce some conjectures about the operators governing 
the finite-size corrections of the Renyi entropies of critical 
quantum chains with PBC. We are going to state them 
separately 
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A M = M = 1/6 M = 1/4 M = 3/10 



0.5 


1.9 


2.0 


2.0 


2.1 


^/2/2 




2.0 


2.0 


2.1 


0.9980 




1.98 


2.0 


2.1 


2.0 




1.97 


2.0 





Table V: The exponent v of the finite-size corrections of the 
von Neumann entropy for the spin-1/2 XXZ chain with PBC, 
and some values of A and magnetization per site M . We 
extracted the exponent, as we did for the other models, con- 
sidering systems with size L = 96. For value of anisotropy 
close to the isotropic antiferromagnetic point (A = 1) and 
M = 0, we were not able to obtain the exponent, mainly due 
to the large finite-size corrections. 

(a) Entropy oscillations in the a-Renyi entropies of 
critical quantum chains only appear for a > 1 and for 
models having at least one [/(I) symmetry, like the spin-s 
Heisenberg, t — J and Hubbard models. Models possess- 
ing only discrete symmetries like the Ising, Blume-Capel 
and Potts models show no oscillations for any value of a. 

(b) For any critical quantum chain the amplitude of 
the leading correction of the a-Renyi entropy with a > 1 
has a universal power decay Pa = 2Xg/a, where is the 
dimension of the energy operator of the model. We stress 
that these results are expected only for the entanglement 
entropy of a single interval. In the case of two disjoint 
intervals, the results of Ref. [s^ indicate that instead of 
the energy operator the leading contributions comes from 
nonlocal spinor operators. 

The amplitudes of the leading corrections for a = 1, 



have a quite distinct behavior from the a > 1 cases. For 
all models we considered these amplitudes have a power 
law decay £^'^ , where ly — pi ^ 2. Since these corrections 
for the a = 1 Renyi entropy (von Neumann entropy) 
are not known in the case of the spin-1/2 XXZ chain 
[anisotropy A] , we present in Table |V] those corrections 
for some values of anisotropy and magnetization per site 
M. We clearly see, that also for the XXZ quantum chain 
the leading finite-size correction always decay as In 
the case of the XX quantum chain, where = 2 is an 
exact resultfi^ or in the XXZ chain we cannot identify 
the operator responsible for such finite-size correction. 
All these results indicate the following conjecture. 

(c) The leading finite-size corrections of the a = 1 
Renyi entropy, or the von Neumann entropy, of any quan- 
tum chain decays as This is the main contribution 
coming from the conical singularities in the conformal 
mapping, but we do not know what operators produce 
such universal correction. 

We expect that the above conjectures can be confirmed 
for other models and we hope they can be understood by 
using the general properties of the underlying CFT ruling 
the long-distance physics of quantum chains. 
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